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Local scale invariance (LSI) has been recently proposed as a possible extension of the dynamical 
scaling in systems at the critical point and during phase ordering. LSI has been applied inter alia 
to provide predictions for the scaling properties of the response function of non-equilibrium critical 
systems in the aging regime following a quench from the high-temperature phase to the critical 
point. These predictions have been confirmed by Monte Carlo simulations and analytical results 
for some specific models, but they are in disagreement with field-theoretical predictions. By means 
of Monte Carlo simulations of the critical two- and three-dimensional Ising model with Glauber 
dynamics, we study the intermediate integrated response, finding deviations from the corresponding 
LSI predictions that are in qualitative agreement with the field-theoretical computations. This result 
casts some doubts on the general applicability of LSI to critical dynamics. 
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I. INTRODUCTION 

The non-equilibrium dynamics of systems at the criti- 
cal point has been recently the subject of a renewed in- 
terest due to its similarities with the aging phenomena 
occurring in a variety of glassy materials, see Refs. Ill2l3l 
for reviews. In the typical scenario the system is pre- 
pared in the high-temperature equilibrium state. Then, 
at time t = 0, it is quenched to the critical tempera- 
ture T c . Due to the critical slowing down, after this sud- 
den thermal perturbation the system undergoes a nev- 
erending relaxation towards the equilibrium state. In- 
teresting scaling properties of this evolution are revealed 
by two-time quantities such as the correlation function 
C x (t,s) = {(f> x {t)(f>o(s)) and the linear response function 
Rx{t,s) = S((f> x (t)) /5h (s)\h=o, where <j) x (t) is the order 
parameter at the point x and time t and ho(s) is the 
conjugated field at time s < t and position x = (here 
translational invariance is assumed). 

In thermal equilibrium C and R depend on t — s 
and they are related by the fluctuation-dissipation the- 
orem, whereas in the present case they are non-trivial 
functions of both times. This is the so-called aging 
regime In particular, in the dynamical scaling regime 
s ^> Tmicro and t—s ^> r m j cro , where T m i cro is some micro- 
scopic time, the autoresponse function i? x =o is expected 
to scale asi£ 

R*=o(t,a) = A R {t-s)-^ +1 \t/s) e T R {s/t). (1) 

The scaling function J-r(v) is universal in the sense that 
it depends only on the universality class of the statistical 
system, provided one fixes the non-universal constant Ar 
to have J-r(0) = 1. The exponent 6 is related to the 
initial-slip exponent of the magnetization 4 (equivalently 
to the autocorrelation exponent^ A, as 6 = 1 + a — X/z), 
whereas a = 2(3 /(vz) — (d— 2 + rf)/z. Here f3, v and r\ 



are the usual static critical exponents, d is the number of 
space dimensions and z is the dynamic critical exponent. 

Eq. (JIJ is basically obtained (for example within 
the renormalization-group approach^) by exploiting the 
fact that, at the critical point, the response function 
transforms covariantly under dilatations in space x t— ► 6x 
and time t i— > b z t. It is possible to prove, under quite 
general assumptions, that the covariance of static cor- 
relation functions (of suitable operators) under a scale 
transformation x i— > fox extends to conformal transfor- 
mations, resulting in strong constraints on the functional 
form of the associated scaling functions. It is then natural 
to address the question whether an analogous extension 
carries over to the case of dynamic scaling with z ^ 1, 
where the dilatation factors in space and time are differ- 
ent. (In strongly anisotropic systems the role of the time 
is formally played by a special direction in spaced.) This 
possibility has been recently investigated, leading to the 
identification of a group of local scale transformations 
which generalize the dynamic scalingA& The covariance 
of the dynamic correlation functions (of suitable opera- 
tors) under such a group is called Local Scale Invariance 
(LSI) . The conditions under which the dynamical scaling 
implies LSI have not yet been established beyond the spe- 
cific case^ z = 2. In the absence of such a general result 
one can heuristically assume LSI and then check whether 
its predictions are actually confirmed. For the response 
function in aging systems [see Eq. JTJ] LSI yieldsZ&iii 

*£ SI >(«)=1. (2) 

This prediction has been found in excellent agree- 
ment with Monte Carlo simulations of two- and three- 
dimensional Ising models and three-dimensional XY 
model, with Glauber dynamics, both at and below 
!T ci 7 i ii i i2 4 i3 In all these simulations the measured quan- 
tity has been either the autoresponse function i? x =o {t, s) 



2 



or some related integrated responses. Moreover the 
prediction @ has been also confirmed by analytical 
results for spherical models with non-conserved order 
parameter^^ and for the two-dimensional XY model in 
the spin- wave approximation. 15 These findings strongly 
support the idea that LSI is actually realized at the crit- 
ical point, providing constrains on the form of scaling 
functions. 

For the spatial dependence of the response function 
LSI predicts^ 



R x (t, s) = R x=0 (t, a)$(|x| (t - S y l ' z ) , 



(3) 



where the function <E> satisfies a (fractional) differential 
equation whose general solution is reported in Ref. |8|. 

Instead of studying numerically local quantities, one 
can consider global ones such as the time-dependent to- 
tal order parameter M(t) = J d d x(j) x (t) = </> q= o(i) (q 
is the wave- vector) . The linear response of M(t) to an 
homogeneous field H(s) applied at time s < t is given 
by 5{M(t))/5H(s)\ H =o = i? q =o(M) = Jd d xR x (t,s). 
This quantity, for non-conserved dynamics, is expected 
to scale asm 



R<!=o(t,s) = A R (t 



- ,s)-( a+1 )+ d / J 



(t/s) 9 F R (s/t), (4) 



where the non-universal constant Ar, is fixed requiring 
•^(0) = 1 an d Fr{v) is a universal scaling function. 
From Eq. JSJ) one can work out the prediction of LSI for 
this susceptibility, finding 



T 



« LSI »(«)e1 



(5) 



and A 



(LSI) 



Ar J d d x $(x). In spite of the support to 
the LSI prediction Eq. @ provided by the result of Monte 
Carlo simulations of some critical lattice models, recent 
field-theoretical (FT) computations&iSi suggest the pres- 
ence of corrections to Eq. @. For the universality class 
of models in d dimensions with iV-component order pa- 
rameter, short-range interactions with O(N) symmetry, 
and purely dissipative dynamics (model A of Ref. Il7f) . it 
has been predicted 



T R FT \v) = l-e 2 c N Af(v)+0(e 3 ) 



(0) 



where c N = 3(N + 2)/[8(N + 8) 2 ], e = 4 - d > 0, and 
A/(«) = /(0) - f(v) > (0 < v < 1) is a monoton- 
ically increasing and regular function of order of unity, 
whose expression can be found in Ref. E In particular 
Eq. © with N = 1 provides, in the dimensional ex- 
pansion, an analytic prediction for the universal scaling 
function of the Ising model with Glauber dynamics. Note 



that J- R r 1> = T R LSI ^ +0(e 2 ), i.e., the prediction of LSI co- 
incides with the Gaussian approximation for the scaling 
function. This observation accounts for the agreement 
between LSI and the analytical results for spherical mod- 
els with 2 < d < 4, mentioned above. Indeed, apart from 
non mean-field exponents, they typically display Gaus- 
sian universal scaling functions. 



The discrepancy between J-r T ^ and J-r SI ^ casts some 
doubts on the effective realization of LSI at the critical 
point and calls for a more careful analysis of the results 
obtained by Monte Carlo simulations, given that no dis- 
crepancy is emerging from the data presented in the past. 



II. THE INTERMEDIATE INTEGRATED 
RESPONSE 



We consider the susceptibility (s < t) 



a/1 



du Rq=o(t, u) 



(7) 



which is the integrated linear response of the order pa- 
rameter to a spatially homogeneous field switched on 
during the time interval [s/2,s]. This integrated re- 
sponse, corresponding to the intermediate protocol pro- 
posed in Ref. 0, presents some advantages over the 
more commonly studied thermoremanent susceptibility 
PTRM(t,s) = fydu Rq=o(t, u) and zero-field cooling sus- 
ceptibility xzFc(t,s) = J^du Rq=o{t, u). As first noted 
in Ref. Il8l ptrm is in general hampered by the presence 
of a finite-time correction which originates from the re- 
sponse to a change in the initial condition. Indeed at the 
lower integration limit the necessary conditions for the 
dynamical scaling regime are not fulfilled and the scaling 
function J2J can not be used in that regime. As discussed 
in Ref. |lj| a similar remark applies to xzfc close to the 
upper integration limit, which contributes with a leading 
term that is independent of the waiting time s, the ex- 
pected scaling part being only a sub-leading correction. 
These problems with the applicability of the scaling form 
(|4"|) are not encountered for the intermediate integrated 
response {7}. 

The expected scaling behavior of x(t, s) can be worked 
out from Eq. |@J: 



(8) 



where A x = ArBq, 
T x {v) = Bg 1 [ dw(l- vw)- {a+ ^ +d ' z w- e T R {vw) , 

Jl/2 

(9) 

and B = (1 - 2 e - 1 )/(l - 9), so that T X (Q) = 1. The 
prediction of LSI that follows from Eq. (0) is given by 

D-l 

^ LSI) (v) =Y i - e [2^1(1 -0,1 + a- d/z, 2-9-v) 

- 2 e - 1 2 F 1 (l -9,1 + a- d/z, 2-9; v/2)] 

(10) 

whereas the FT prediction is obtained from Eq. iJBJ, with 
9 = O(e), z = 2 + 0(e 2 ), and rj = 0(e 2 ): 



= l-e 2 c N - [ duAf(u) + 0(e 3 ) (11) 

v J v /2 
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We study the intermediate integrated response J7J by 
simulating Ising models on square and on cubic lattices. 
The systems are thereby prepared in a completely uncor- 
rected initial state and then quenched at time t = to 
the critical point. The temporal evolution is realized us- 
ing the standard heat-bath algorithm. At time t — s/2 a 
spatially constant field with strength H = 0.05 J (J being 
the strength of the nearest neighbor couplings) is applied. 
This external field is switched off at t = s. We checked 
that for this field strength we are well within the linear 
response regime, as shown in the inset of Figure QJi. The 
data discussed in the following are free from any finite- 
size effects and have been obtained for systems consisting 
of 450 x 450 spins in two dimensions and 80 x 80 x 80 spins 
in three dimensions. Typically, we averaged over 50000 
different runs with different realizations of the thermal 
noise. 

The temporal evolution of the integrated response 



1 



i=l 



t > S 



(12) 



is shown in the insets of Figure 1 in two and three di- 
mensions. Here Si(t) G { — 1, +1} is the value of the Ising 
spin located at the lattice site i at time t. The sum in 
(|12|l is over all n sites and (...) indicates the average over 
the thermal noise. A similar behavior is observed in two 
and three dimensions: for a fixed value of the waiting 
time s the susceptibility rapidly displays a power-law in- 
crease as a function of t/s, with exponent (d — X)/z [see 
Eq. |JSJ]. Its measured values in two and three dimen- 
sions (0.186(2) and 0.108(2), respectively) are in excel- 
lent agreement with the expected values 0.185 and 0.108, 
obtained by using the known values of A (d = 2: 1.60, 
d = 3: 2.78) and z (d = 2: 2.17, d = 3: 2.04)££ 

The scaling behavior JHJ is tested in the main images 
of Figure 1. Plotting t a ~ d ' z x as a function of s/t yields 
a remarkable data collapse. The scaling functions so ob- 
tained can be compared with the available analytical pre- 
dictions. In particular, assuming the well-known values 
of the critical exponents appearing in Eq. (JSJ one de- 
termines the non-universal constant A x from the small 
s/t behavior of t a ~ d / z \- Then the prediction of LSI and 



e£{LSI),{FT) 



where T- 



(LSI) 



is 



FT are given by A x {s /t) 1 " 9 T] 
provided in Eq. (JTOJl with 6 = 0.38, a = 0.115 in d = 2 
and 9 = 0.14, a — 0.506 in d — 3. The numerical estimate 
of F X FT) is obtained through Eq. JTTJ with N = 1, e = 2 
and 1 in two and three dimensions, respectively, corre- 
sponding to its Pade approximant [2,0], which does not 
differ significantly from the approximant [0, 2]. The com- 
parison shows that the LSI curves (black lines) system- 
atically lie above the numerical data for larger values of 
s/t. Including the two- loop correction coming from field 
theory shifts the theoretical curves (grey lines) closer to 
the data. It has to be stressed that this is the first time 
that the existence of corrections to LSI in critical non- 
equilibrium systems has been numerically proven. 




FIG. 1: Rescaled intermediate integrated response vs s/t for 
different waiting times s computed in the critical Ising model 
in (a) two and (b) three dimensions. In both cases a perfect 
data collapse is observed. The LSI prediction lies systemati- 
cally above the numerical data for larger values of s/t, whereas 
the two-loop correction shifts the theoretical curve closer to 
the data. Error bars are smaller than the sizes of the symbols. 
The insets display the power law increase of \ as a function 
of t/s for fixed values of s. In the inset of (a) we also probe 
the linear response regime. Solid lines have been obtained 
with the field strength H — 0.05J, whereas the symbols for 
s = 200 and s = 100 result from simulations with H = 0.025J 
and H = 0.1J, respectively. The measured susceptibility does 
not depend on H, showing that we are well within the linear 
response regime. 



Finally, we show in Figures 2a and 2b the scaling func- 
tion jF x (s/t) obtained from the data of Figure 1 after a 
multiplication by (s/t) 0-1 and a proper normalization. 
Corrections to the LSI predictions are clearly revealed in 
this quantity both in two and three dimensions. Analyz- 
ing the scaling function T x (s /t) seems therefore to be the 
appropriate way for highlighting differences between the 
theoretical predictions and the numerical data. In the 
more common approach where the rescaled susceptibil- 
ity s a - d / z x is plotted against t/s corrections are hardly 
detectable in three dimensions. 
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FIG. 2: Comparison of the numerically obtained scaling func- 
tions T x with the theoretical predictions in (a) two and (b) 
three dimensions. Corrections to LSI are clearly observed for 
this quantity in both cases. In (a) error bars are comparable 
to the symbol sizes. 

III. CONCLUSIONS 

Local Scale Invariance proposes to generalize dynami- 
cal scaling to local space- and time-dependent scale trans- 



formations, leading to several predictions for dynamical 
correlation and response functions. Some of these pre- 
dictions have been verified in the past through exact so- 
lutions of soluble models and through numerical simula- 
tions of non-soluble ones. In this paper we have focused 
on the scaling functions of the dynamic response function 
in non-equilibrium critical Ising models. Previous nu- 
merical studies did not reveal any systematic deviations 
from the LSI prediction when computing local quantities. 
Recent field-theoretical calculations, however, yielded a 
two-loop correction to the LSI prediction. In the present 
work we have studied the non-equilibrium response of 
the total order parameter to an homogeneous external 
field, i.e., we have investigated a global quantity. It has 
been realized recently that in the usually studied inte- 
grated responses, the thermoremanent susceptibility and 
the zero-field cooling susceptibility, additional terms ap- 
pear that make the study of the dynamical scaling part 
(i.e., the aging part) notoriously difficult. We therefore 
propose to study the intermediate integrated response 
where the field is only switched on during a time inter- 
val [s/2,s], with < s < t, t being the time at which 
the resulting magnetization is measured, and the time 
at which the quench to the critical point occurs. Look- 
ing at this quantity we have indeed identified corrections 
to the LSI prediction. In addition, we observe that the 
field-theoretical two-loop correction brings the theoreti- 
cal curve closer to the numerical data. 
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